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1 Introduction 



The Gribov-Zwanziger framework [H [21 E] consists in restricting the domain of integration in the Feyn- 
man path integral within the Gribov horizon. It has motivated extended studies of nonperturbative tools 
for investigating the infrared behavior of the gluon and ghost correlation functions, see [I] for the use of 
modified Schwinger-Dyson equations for QCD, El \7\ [8] for recent analytical results and [91 [TUl EE] for 
numerical data obtained through lattice simulations. 

Zwanziger has been able to show that the restriction of the path integral within the Gribov horizon 
for the gluon can be achieved by adding to the Fadeev-Popov action in the Landau gauge a local action, 
depending on new fields, with well-defined interactions with the gluons and their Faddeev-Popov ghosts 
in the Landau gauge [21 [3]. However, this local action violates the BRST symmetry by a soft term, 
so that BRST symmetry is only enforced in the scaling limit. The yet unorthodox point of view that 
the BRST symmetry can be broken in the IR region of QCD was heuristically anticipated by Fujikawa 
|12j . It is not in contradiction with any given physical principle, since it is by no means necessary that 
the QCD microscopic theory possesses a unitary sector for its partons, namely the quarks and gluons, 
to warrantee unitarity properties among the sector of bound states that constitute its spectrum. The 
only requirement is that the modified theory remains renormalizable and that the BRST symmetry is 
recovered in the ultraviolet region, in order to suitably describe the asymptotic properties in terms of 
almost deconfined partons, as predicted by asymptotic freedom and short distance expansion. From a 
physical point of view, it is gratifying for the quark-gluon model introduced in [13] that the modification 
of the usual Feynman propagators into a Gribov-type propagators eliminates from the beginning all par- 
tons from the spectrum, since their modified propagators have no poles on the real axis, a property that 
anticipates quite well the confinement. 

What justified our previous work [13] is that the genuine geometrical approach of Zwanziger leaves 
aside the quarks, which do not participate to the Gribov phenomenon, while the idea of a parton model 
suggests that quarks and gluons should be treated on the same footing, and quark propagators should 
have an analogous behavior as the gluon one in the infra-red domain. An idea for getting such modified 
quark propagators was thus needed, which goes beyond the Gribov question. A generalization of the 
work of Zwanziger was also needed [61 [8] in order to achieve a different behavior for the shape of 
the genuine Gribov-gluon propagator D(q 2 ) that vanishes at q 2 = 0, for accommodating recent lattice 
simulations that seem to indicate that the gluon propagator goes to a non-vanishing constant at very 
small q 2 PEIilCLl]. 

To provide such an improved, local and renormalizable, quantum field theory that gives the wanted mod- 
ifications in the infra-red region, both for quarks and gluons, we proposed in [13] the following picture. 
Given a theory of partons (eg quark and gluons in four dimensions, or a scalar fields in dimensions such 
that the ultraviolet divergences are renormalizable or super-renormalizable), it can always be coupled to 
a topological field theory made of new fields arranged as BRST trivial doublets, in such a way that the 
partons are already confined at the tree level by their mixing with the unphysical fields of the topological 
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field theory. We found that such a mixing can be generally allowed by a soft breaking of the BRST 
symmetry. Here, confinement is meant in a very simple way. The propagators of all fields have only poles 
at complex positions^- This implies that the theory has no vacuum for the partons and all observables are 
made of composite operators, defined by solving the cohomology of the BRST operator. These composite 
operators can be renormalized in the standard way, with expectation values related to the parameters of 
the soft breaking mechanism. This idea was inspired by the algebraic characterization of the local terms 
that Zwanziger introduced to complete the Faddeev-Popov action for the gluons. Part of the Zwanziger 
action can be recognized as a topological action, which involves bosons and fermions that transform under 
the BRST symmetry as a system of two BRST trivial doublets, and it is BRST-exact. The remaining 
part of Zwanziger's action breaks the BRST symmetry in a soft way, and yields a Gribov propagator for 
the transverse gluon. 

The addition of the new fields arranged as BRST-exact doublets eventually allows for the introduction 
of massive parameters that can be used to modify the infrared behavior of the theory, without changing 
the set of observables. The necessity of breaking the BRST symmetry can be easily understood within 
this framework. If the added action were BRST-exact, nothing would be changed for the predictions of 
the original parton theory. Indeed, the observables, defined as the cohomology of the BRST operator, are 
the same and the effect of integrating over the new fields would only consist in multiplying the partition 
function by one. So, if the new field dependence is through BRST-invariant terms (and thus BRST-exact 
terms because they transform as BRST-exact doublets), there is no way to improve the infra-red behavior 
of the amplitudes. One can notice that the possibility of adding such a BRST-exact terms can be related 
to the formal invariance of the path integral under arbitrary changes of variables for the parton fields. 
By doing such changes of variables and playing with Lagrange multipliers and determinants formula, 
one can indeed recover the class of BRST-exact action that are bilinear in the new fields that have the 
opposite statistics to the partons^). To obtain the wanted modification of the infra-red behavior of parton 
correlators, one must go further, and consider the possibility of an explicit soft breaking of the BRST 
symmetry. This, of course, can modify the parton propagators in the low energy domain, without chang- 
ing their ultraviolet behavior, provided the theory is renormalizable. This will also modify the numerical 
values of the observables, by giving them a dependence from the parameters of the soft breaking. 

In the case of the improved quark-gluon theory, the properties of the breaking and the respective conse- 
quences for the theory can be summarized as follows: 

• The breaking term is soft, meaning that its dimension is smaller than the space-time dimension. 
As a consequence, the breaking can be neglected in the ultraviolet region, where one recovers the 
notion of exact BRST symmetry. 

• The BRST operator preserves nilpotency. Moreover, the set of physical operators of the theory, 
identified with the cohomology classes of the BRST operator, is left unmodified. This occurs 

1 More precisely, the propagators display violation of reflection positivity, a feature which invalidates the interpretation 
of partons as excitations of the physical spectrum of the theory. 

2 In string theory, this idea has been already used by adding so-called "topological packages" to any given 2d-string action, 
which allows one to show in this way possible relationships between the different string models [141 1151 HE] . 
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because the new additional fields are introduced as BRST doublets. As such, they do not alter the 
cohomology of the BRST operator. 



• The soft breaking of the BRST symmetry is introduced in a way compatible with the renormal- 
izability. More precisely, it is associated to a quadratic term in the fields obtained by demanding 
that the additional fields couple linearly to the original fields of the theory, so that the resulting 
propagators turn out to be modified in the infrared region. As a consequence, a given correlation 
function can display a different behavior, when going from the deep ultraviolet to the infrared region 
in momentum space. 

• The soft BRST symmetry breaking is meant to be an explicit breaking, i.e. it is not a spontaneous 
symmetry breaking, which would give rise to Goldstone massless fields, and thus to a quite different 
framework, as it will be explained elsewhere. 

• All massive parameters should be related to the unique scale of the theory, namely Aqcd ; by 
requiring that the massive soft parameters satisfy suitably gap type equations which allow to 
determine them in a self-consistent way. This is the case, for example, of the massive Gribov 
parameter 7 which is fixed by a gap equation [21 El El El E] , see Sect. 4. 

Following the procedure described above, in [13] a model accounting for a soft BRST symmetry breaking 
giving rise to a modification of the long distance behavior of the quark propagator was established. In 
fact, the model predicts 



is a function depending on the soft breaking mass parameters M\, M2 and m. From expression ([2]) one 
sees that the function A(k) vanishes in the deep ultraviolet region where the usual perturbative behavior 
for the propagator is recovered. It is worth mentioning that expression ([2]) provides a good fit for the 
dynamical mass generation for quarks in the infrared region in the Landau gauge, as reported by lattice 
numerical simulations of the quark propagator [171 I18j . We recall that the function ([2]) is analogous 
to the one appearing in the gluon propagator within the Gribov-Zwanziger framework [21 E]- In fact, 
as discussed in |13j . the Gribov-Zwanziger action can be recovered through the introduction of a soft 
breaking, related to the appearance of the Gribov parameter 7. This parameter is needed to implement 
the restriction to the Gribov horizon, which turns out to be at the origin of the soft breaking of the 
BRST symmetry [HIE]- In particular, for the tree level gluon propagator, one finds 




(1) 



instead of the standard Dirac propagator for the quarks ^ M . In this equation 




(3) 



with 



M 2 (k) 



k 2 + [i 



2 ' 



(4) 
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where the second mass parameter \i accounts for the nontrivial dynamics of the auxiliary fields needed 
to localize Zwanziger's horizon function, see [HE]. 

It is of course striking how similar are the modified gluon and quark propagators, in the sense that 
the introduction of the massive parameters by the above mechanism has eliminated the infra-red prob- 
lem for both type of fields. Within the context of plain perturbation theory, these parameters can be 
fined-tuned, and even be used as plain infra-red regulators that are controlled by softly broken Ward 
identities. However, in a confining theory like QCD, they must be related to physical observables, and 
non-perturbatively computed as a function of the basic parameter of the theory, Aqcd- The clear ad- 
vantage of our approach is that the modified theory is no more submitted to the Gribov ambiguity and 
the parton (quark and gluons) propagators have no poles at real values. The model excludes their ap- 
pearance in the spectrum from the beginning, and its perturbative expansion can be possibly compared 
to the prediction of non-perturbative approach such as the lattice formulation in the Landau gauge. The 
method solves in a quite simple conceptual way the necessity of introducing massive parameters, not only 
for the Gribov ambiguities but also for defining the scale of the chiral symmetry breaking. 

The goal of this paper is to prove an elementary, but necessary property of this model: it has to be 
multiplicatively renormalizable, to ensure that the breaking of the BRST symmetry remains soft, at any 
given finite order of perturbation theory, and that no further parameters than those which modify the 
infra-red behavior of propagators can appear. We will do this by employing the algebraic renormalization 
|19j . Notice that such a proof can be generalized to other cases, where one may wish to introduce an 
infra-red cut-off for other theories with nontrivial gauge invariance and control its effect by softly broken 
Ward identities, within the context of locality. This could be of relevance for the study of non-exactly 
solvable super-renormalizable theories, like 3-dimensional gauge theories, where one must control the 
way the infra-red regularization is compatible with the BRST symmetry, or even for gravity, when it is 
computed at a finite order of perturbation theory. 

The paper is organized as follows. In Section [2] we provide a short overview of the model and of its 
soft BRST symmetry breaking. In Section [3] the algebraic proof of the renormalizability is given. In 
Section [4] we consider the inclusion of the Gribov- Zwanziger action and we discuss the renormalizability 
of the resulting model. Finally, in Section [5] we collect our conclusions. 

2 The quark model and its BRST soft breaking term 

We start with the Yang-Mills action quantized in the Landau gauge, 




) 



(5) 



where 




(6) 
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is the covariant derivative in the fundamental representation of the SU(N) gauge group, with generators 
{T a yi, and 

D« b = d ll 5 ab -gf abc A^, (7) 

is the covariant derivative in the adjoint representation. The first set of small Latin indices {i, j, . . .} E 
{1,..,N} will be used to denote the fundamental representation, while the second set {a, b, ...,h} € 
{1, .., N 2 — 1} will be employed for the adjoint representation. The set of Greek indices {a, (3, 7, 5} stand 
for spinor indices. The remaining Greek indices will denote space-time indices. 

The action ([5]) is invariant under the BRST transformations, namely 





= -D ab c b 






= -igc a {T a ) ij 






= -ig^ a c a (T a 




sc a 


= ^gf abc c b c c , 






= ib a , 




sb a 


= . 





Following [13], a model for the dynamical quark mass generation can be constructed by introducing 
two BRST doublets of spinor fields (£*, 9 l ) and (rf, A*), transforming as 

s& = & a , < = 0, 

sr f a = AL, sX^ = . (9) 
The propagation of these fields is described by the following BRST-exact action 

% = s f d A x{-ff a d 2 e a + e a d 2 ri a + m\ff a e a -eM) 

= j d*x (-Kd 2 & - e a d 2 K - f? a d 2 V a + W a d 2 rf a + m 2 (A*X + C a \\ + ffji - «)) , (10) 

where m is a mass parameter. Further, we introduce the coupling of the spinors and {rf , A*) with 

the matter field ip l a , 

S M = j d 4 x (Mj^M + - M 2 (A^< + fcX'j) . (11) 

Evidently, the action (|10p is BRST-invariant. This is not the case of Sm, which will give rise to a soft 
breaking of the BRST symmetry, parameterized by the two soft mass parameters M\, In fact, 

sS M = A , (12) 
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where 

+ igM 2 c a (Xl t (T a ) ij ^ a - fi a (T a f\i)) . (13) 
Notice that A is a soft breaking, i.e. it is of dimension less than four. 

Performing the integration over the auxiliary spinor fields, yields a nonlocal action for ip and 

S f = J d*x ( 7M ) a/3 DjtyJ - 2M?M 2 ^ a (^^2) C) • (14) 

It is important to emphasize that the introduction of the auxiliary spinor fields can be seen as a tool for 
localizing the nonlocal term (|14j) that appears in the fermionic sector, in the same way as the nonlocal 
Gribov-Zwanziger horizon function is cast in local form through the introduction of a suitable set of 
auxiliary BRST doublet fields [21 EJ El [7] . 

2.1 Introducing sources for controlling the soft symmetry breaking 

In order to prove that the model described by the action 

S = S mv + S^x + Sm , (15) 

is multiplicatively renormalizable, we follow the procedure outlined by Zwanziger in the study of the 
Gribov horizon in the Landau gauge [21 [3]. It amounts to embedding the action (|15p in a larger model, 
S — > So, displaying exact BRST invariance. This is achieved by treating the breaking term A as a 
composite operator, which is introduced into the theory through a suitable set of external sources. The 
original action S is thus recovered from the extended action So by demanding that the sources acquire a 
particular value, which we shall refer to as the physical value. The renormalizability of S follows thus by 
proving the renormalizability of the extended action Sq. 

In order to introduce the extended invariant action So, we make use of the following set of sources 
(J,H), (J, 8), (K, G), (K, G) and (N,P), assembled in BRST doubletfl i.e. 





- H lJ 


brL a/3 


= 




- H ij 


bn aj3 


= 






SLr a/3 


= 0, 


SJX a(3 


- U af3 > 


SU a/3 


= 0, 


sN 


= P, 


sP 


= 0. 



(16) 



3 For future purpose, we recall that all quantum numbers of fields and sources are displayed in Tables [T] and [2] including 
the charge Q4N that will be defined through expression (123[) . 
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The invariant action that accounts for the extra spinor fields and for the breaking term is thus defined as 

SjK = 8 J d A X {-ff a d 2 ti + + P(ff a e a - CM + *NP 

which, explicitly, reads 

Sjk = J d A x (-A^ 2 Ci - C a d 2 X a - ff a d 2 6\ + & a d 2 rf a + P(\Ui + £ Q A Q + ffji - «) + aP 2 

+ G%\^ + G%^\\ + K%X a igc* (Ty k - K%ig^c a {I*)* A*,) , (18) 

whose BRST invariance is manifest. The parameter a, in expression (|18p . is a dimensionless parameter 
needed for renormalization purposes. For the so called physical values of the external sources, we have 





phys 


— J a/3 


phys 


= o, 


K ij 

■"■a/3 


phys 


~ ""a/3 


phys 


= o, 


TTij 

n a/3 


phys 




phys 


= M 2 S ij d a /3 , 


^a/J 


phys 


~ ^a/J 


phys 


= -M 2 S*6 afi 


P 


phys 


= m 2 , 






N 


phys 


= . 







(19) 



The BRST-invariant extended action So is thus defined as 

So = Sinv + + Sjk + S ext . (20) 

It is easily checked that the starting action S, eq. (|15|) . is recovered from the extended action Sq when 
taking the physical values of the sources, eq.(|19|). namely 

= S + [ d A x am 4 . (21) 
phys J 

It is worth noticing that the action So possesses an additional U (47V) global symmetry, provided by 

Q%S = , (22) 

where 

n ij f , 4 / £j 5 rj S j 5 ^rj 5 j 5 ai ^ j $ -j $ jjk $ 

Qap ~ j r^r w P p s%~ p m + a w p + ^w a ~^w a + Jf *8JE 

jjk 5 ik 5 _ - ik 5 ik 5 _ - ik 5 jk 5 _ - ik 5 \ 
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fields 


A 


b 


c 


c 


iji 


ijj 


£ 




A 


A 


9 


9 




i] 


UV dimension 


1 


2 





2 


3/2 


3/2 


1/2 


1/2 


3/2 


3/2 


1/2 


1/2 


3/2 


3/2 


Ghost number 








1 


-1 




















1 


1 


-1 


-1 


Q4N - charge 




















1 


-1 


1 


-1 


1 


-1 


1 


-1 


Spinor number 














1 


-1 


























Statistics 








1 


-1 


1 


-1 


1 


-1 


1 


-1 


2 








-2 



Table 1: Quantum numbers of the fields. 



sources 


n 


L 


Y 


Y 


J 


J 


K 


K 


H 




G 


G 


N 


P 


UV dimension 


3 


4 


5/2 


5/2 


2 


2 


1 


1 


2 


2 


1 


1 


2 


2 


Ghost number 


-1 


-2 


-1 


-1 


-1 


-1 


-1 


-1 














-1 





QiN - charge 














1 


-1 


1 


-1 


1 


-1 


1 


-1 








Spinor number 








1 


-1 


1 


-1 


1 


-1 


1 


-1 


1 


-1 








Statistics 


-1 


-2 





-2 


1 


-3 


1 


-3 


2 


-2 


2 


-2 


-1 






Table 2: Quantum numbers of the sources. 



The trace of the operator Q*^g, i.e. Qf" = Q$n, defines a new conserved charge, and thus an additional 
quantum number for the fields and sources, allowing for the introduction of a multi-index I = (i, a), with 
I £ {1, ... , 4iV}, for the fields 4> = (£> V> K 8) an d sources V = (K, J, H, G). Accordingly, we shall set 

C = ¥ , 

1$ = r ?'" (24) 
From now on, we shall make use of the multi-index notation. 



3 Algebraic proof of the renormalizability 

Let us face now the issue of the renormalizability of the extended action So, a task which we shall 
undertake by making use of the algebraic renormalization [19]. Let us start by establishing the set of 
Ward identities fulfilled by the action So. To that purpose, we add external sources, L a , $7^, Y^ and Y£, 
coupled to the non- linear BRST variations of the fields, namely 

s ex t = s j 'd 4 x(^ + iv + yx + ^) 

= J d A x (-n^D^ + l - 9 r hc L a c b c c - igYicf (T a f < - ig¥ a c a (T a f Y<\ , (25) 

with 

stt^ = sL a = sY = sY = . (26) 

The term (|25p allows one to convert the BRST symmetry into the corresponding Slavnov- Taylor identity. 
Thus, to prove the renormalizability of the model we shall consider the more general action 

£ = Sq + S ex t . (27) 
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3.1 Ward Identities 



The complete action (|27p fulfills a rich set of Ward Identities, namely: 
• The Slavnov- Taylor identity 

5(E) = , 

where 

, 4 / 5E <5E 5E 5E 5E 5E <5E i5S 5E - 7 5E 



5(E) y d x ^ ^ ^ + ^ ^ ^ ^ + §ca gLa + i6 a § _ a + tf _ 7 + A 

+ 01 ^ + A ^ + + G ^ + + ^'^J + P §) 

St 1 fiv 1 5Ji J 5Ki 3 5JL J 5Ki J 8N J 



• The gauge condition: 

• The antighost equation: 

• The ghost equation: 
with 



— - ir) A a 



<5E SE 

6c* + "sm u 



S a S = AS 



and 



The classical rigid invariance 



where 



# 7J E = , 



>U / ,j4„. / \7 ° \J ° , rj-Ifc ° rjJk fi , 7 Jk fi , jlk 



R ~ d X \ X JU~ X JH +H « JnJk - H « JnTk + J « TWIk + J * 



5£ J 8? a 8G J a k a 8G Ik a 8K lk a 5K J a k 
• The equations of motion of the doublet fields: 



Sfj 1 
8E_ 
W 



5X 1 a 5Yi 
8? + 8Yl 
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The Ward identity for the source N: 



3.2 The invariant counterterm 



SE , , 



In order to characterize the most general invariant counterterm which can be freely added to all orders 
in perturbation theory [19], we perturb the classical action £ by adding an integrated local polynomial 
Yfount Q £ dimension bounded by four, and with vanishing ghost number as well as Qhn charge. We 
demand thus the perturbed action, (E + e $] coun *), where e is an expansion parameter, fulfills, to the first 
order in e, the same Ward identities fulfilled by the classical action E, i.e. eqs. (|28f l- (|38|) . This requirement 
gives rise to the following constraints for the counterterm Y l count : 





= o, 


(39) 


6b a 


= o, 


(40) 


^ i Ft 1 sr-count 


= o, 


(41) 




= o, 


(42) 




= o, 


(43) 


8N 


= o, 


(44) 


^Y coun t 
Srj 1 


= o, 


(45) 


^Y coun t 
S9 1 


= o, 


(46) 




= o, 


(47) 


^Y coun t £Y coun t 
8? ' J ° 5Yi 


= . 


(48) 



where the operator By, in eq. (!39p stands for the nilpotent linearized Slavnov- Taylor operator, 



f 4 / SE S . 5E 5 SE 5 SE 5 SE 5 t 6E 5 

5E S 5E 5 ... 5 m 5 - f 5 nI 5 xJ 6 
5c a 5L a 5L a 5c a <5c a S^ 1 Sfj 1 d^ 1 Srj 1 

6 J* 3 SK* 3 a 5 J 1 3 a 6K* 3 SNJ y ' 



The first constraint, eq. (|39p . identifies the invariant counterterm as the solution of the cohomology of the 
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operator By, in the space of the integrated local field polynomials of dimension four. From the general 
results on the cohomology of Yang- Mills theories [)"!) . it follows that Y, count can be written as 



^count = ^0j d ± xF a vF a v + ^(-1) ? (50) 

where A*-" 1 ) is a local integrated polynomial in all fields and sources, with dimension four, ghost number 
minus one and vanishing Qan charge, namely 

= J (fx (cuA^n* + a 2 d^ a A1 + a 3 c a L a + a 4 (?b a + a^f abc c a c b (f + a % \ l K l J^l 

+ auYiGtjt 1 + ax^J^ii + a^JLH 1 + a^Y* + a 17 Y^ Q + a^S'J j'J + a^PjH^ 

+ am&^K^&i + a 21 c a c a iV + a 22 A a ^N + a^A^J K J J + a 24 G#d 2 ^' + a 2 ^G I jG I jN 

+ a 27 fj T d 2 e + a^dW + a^A^e + a^A^rj 1 + a 31 c a c<W 

+ and^tW + araiW + a^fY + a^e-eR'jG'i + a 36 G T J d 2 + a^PiV 

+ a 38 f A J iV + a 39 AYiV + a^f^N + anff&N + a 42 ^^G^K^) . (51) 

The coefficients {ao, ai, a 2 , • • • , 0,42} in expressions (|50p and (|5ip stand for arbitrary constants parameters. 



After a straightforward analysis using the conditions f)39f )-(|48f) it turns out that the only non-vanishing 
coefficients are: ao, ax, a 2 , ae, a 7 , a 14, ais, aig, 017, 037, with the following relations between them 



and 



Then, after redefining 



a 2 = ax , (52) 



ai6 = axB = a-j , 
-an = ai4 = a 6 . (53) 



a6 + a 7 i-> a 2 , 

a 37 ' ^ cra 3 , (54) 



for the form of the final allowed counterterm one finds 



Z count = J <fx (^*±^±\ F ; v F; v - a 1 d fl c a ^ - ai d,c a d,c a - ai ig^ (l^T^^ 

- a^D^i^^^ + asaP 2 ) , (55) 

which corresponds to the usual Yang-Mills counterterm in the Landau gauge with the addition of an 
energy vacuum term, aP 2 , related to the mass m. 
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3.3 Stability 



It remains now to discuss the stability of the model, i.e. to check that the counterterm Yi count can be 
reabsorbed in the classical action S by means of a multiplicative redefinition of the coupling constant g, 
of the parameter a, anf of the fields and sources [19], namely 

E(g, a, 0, $) + eE count = E(g , a , 0o, *o) + 0(e 2 ) , (56) 

where <j> stands for all fields and $ for the sources, 

4> G {-4, ip, c, c, b, £, £, A, A, 9, 9, 77, fj} , 

$ G {Q,L,Y,Y,J,J,K,K,G,G,H,H,N,P}. (57) 



Thus, by defining 



we obtain 



Z l f = \ + -{a Q + 2 ai ) , 

<7 1 °° 

^ = l-e^, 

V = Z / =1+ 2 a2 ' 
Z CT = l + ea 3 . (59) 

Expressions (|59p constitute the independent renormalization factors. All the remaining factors can be 
expressed in terms of the renormalization factors appearing in eq. (|59j) . In fact, for the Lagrange multiplier 
and Faddeev-Popov ghost fields we have 

— & a 1 

Z\l* = Z]' 2 = Z^Z- 1 ", (60) 
while the renormalization of the external BRST sources are found 

Zn = Z-^Z-'I\ 



^0 = Z l J 2 <j) , 

<3?o = Z$$ , 

9o = Z g g , 

cr = Z a a , (58) 



e 



9 A 

Zl = Z 1 /, 

Zy = Zy = Z^ZTZ- 1 ' 2 . (61) 

As expected from (|55|) . the renormalization properties of the usual Yang- Mills sector are preserved. For 
the doublet fields we obtain 

v l/2 _ ryl/2 _ „l/2 _ ryl/2 _ 

Z f - Zj^ - /j x - A x - 1 , 

Z]' 2 = Zf = Z~ 1 ' 2 = Z," 1/2 = Z\l*ZrW. (62) 
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Filially, for the remaining sources we have 



Z H = Z R = Z G = Z a = Z^ 1/2 , 

Zj = Zj = Z K = Z R = z~ x i 2 z x J^z^ 12 , 

Z P = 1 . (63) 



This ends the proof of the multiplicative renormalizability of the model proposed in this article. For 
completeness, let us give the expression of the bare action written in terms of the renormalized fields and 
parameters: 

S = J d A x(±Z A [d^Al - d v Al) dpA* + Z g zf 9 r mn d,A a v A™A n v + \z 2 g Z 2 A g 2 f abc f amn A\A C V A™ A n v 

+ Z^i ( 7M ) a/3 drfi - iZ g Z^z)l 2 g^ a ( 7/1 ) a/3 A; (T a f ^ + ib a d^ + igY*c a (T a T 3 < 

- ig^ a c a {T a f Yl - Z g l Z A 1/2 (d^ + nj) d^c a + (d^c a + gf abc A c / + ±gL a f abc c b c c 

- A J <9 2 £ Z - £ I d 2 \ I - fj I d 2 e I + J dV + mA J £ J + A J + rmffl 1 - m#V 

+ Mi^lC + Ml$& - M 2 X^i - M 2 CAL + Z a am 2 ^j . (64) 

From expression (|64h one can obtain the renormalized version of (114p . namely 

= J d 4 x (z^i { ltl ) ap - iZ g Z^ 2 g^(^) aP A^T a f^ - 2Z^M 2 M 2 ^ a ^) • 

(65) 



4 Inclusion of the Gribov-Zwanziger term 

4.1 A brief overview of the Gribov-Zwanziger action and of its soft BRST breaking 
term 

The Gribov-Zwanziger framework [TJ EU3] enables one to take into account the existence of the Gribov 
copies, which affect the Landau gaugqj. This is done by restricting the domain of integration in the 
Feynman path integral to the so called Gribov region SI, defined as the set of fields fulfilling the Landau 
gauge condition and for which the Faddeev-Popov operator, Ai ab = —d^D^ b (A), is strictly positive 

SI = {A» , d^Al = . M ah > 0} , (66) 

As shown in [21 E] , the implementation of the restriction to the region SI is done by adding to the starting 
action a nonlocal term, known as the horizon function, namely 

Sgz = -5 V J d 4 x / a6c 4 [(8 ■ Dy 1 ] ad f dec A; . (67) 

4 See [2D] for an introduction to the subject of the Gribov copies. 
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The parameter 7 has the dimension of a mass and is known as the Gribov parameter. It is not a free 
parameter, being determined in a self consistent way through the gap equation 



5T 

$1 



.,-0, (68) 



where V stands for the effective action evaluated in the presence of the horizon function (|67p . Despite of 
its nonlocal character, the term Sqz can be cast in local form by introducing a suitable set of auxiliary 
fields 



e~ bGZ = J DipDif DujDlo e g z , (69) 

+ 57 2 J d 4 x - (pf)Al ~- g (N 2 - 1)7 2 ) • (70) 

Here, (cpt , (pt ) form a pair of complex commuting fields, while {00^, Q^) form a pair of complex anti- 
commuting fields. These fields are assembled in BRST doublets 

.ab , ,ab „, ,ab 



8<f%=U%, «< = 0, 



sQ? = <ptf, spf = 0, (71) 



and, as pointed out in (5j EJ El [13], the local action (j70l) gives rise to a soft breaking of the BRST 
symmetry, due to the presence of the Gribov parameter 7. In fact, it turns out that expression (JTOj) can 
be written as 

ShT 1 = V + 57 2 / d 4 x (V fec « fe - (pf)Al ~- g (N 2 - 1)7 2 ) , (72) 

with 

= f d A x (-<ffd v D?<p* + Q?d v D?u* + (d^gr^^D^^) , (73) 



so that 



„ oLocal .2 a 



A 7 = J d i x[gf^ufA^-gf^f-(pf)Dfc d ). (74) 



In order to keep control of the soft BRST breaking term, we proceed as before and introduce a set of 

jab fjab T/afe t/ a ^ 



external sources (U?£, [/—, V?™, V^) transforming as 



<tV ab = U ab sTJ ab = 



sU$ = V$, sV$ = 0, (75) 
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and whose physical values are defined by 



yab 



phys 
jjab 



_yab 



phys 



phys 
jjab 



phys 



Thus, we can replace the breaking term in eq. (|72j) by the following BRST-invariant source term 



d x U D ip v + V D u u - U V 



ab r\ac,-,cb 



jab t rab 



d 4 x (v#D?<p? - Uf v [B™u? + gf acd (Dffc e )^} 



i rrab r^ac-cb , irab\ nac ,—cb i _ facd / r~ide „e\-cb~\ iTTabjjab -\rab\r 



(76) 



(77) 



Notice that the source term 5 source gives back the original BRST soft breaking term when the sources 
attain their physical values (j76f) . In fact, after a little algebra, one finds 



So, 



phys 



gi 2 f d 4 x (r fec (^-^ fe )^ 



'--(n 2 - i)r 



(78) 



We are now ready to discuss the inclusion of the Gribov-Zwanziger term (|70p into our starting action X, 
eq. (j27|) . To that purpose we consider the more general action 



St 



where 



5/, 



/i 2 s 



^4 _ - ab , ab 



ab ab -ab ab 



(79) 



(80) 



As discussed in [6l [7] , the term takes into account the nontrivial dynamics of the auxiliary localizing 
fields (</?^V^ 6 )- The introduction of the BRST-invariant term follows from the observation that the 
dimension two condensate {(p't <^f 
namely 



. -,ab, ,ab 
Lu,, Lu 



I ,-rflb , ab 



„ / has a nonzero value for non- vanishing Gribov parameter 7, 
2 



-ab ab\ 



m - 1 ) 2 v2 gjV i/2 2 

64tt 



(81) 



The existence of this condensate is taken into account through the mass parameter \i which, in a way 
similar to the Gribov parameter 7, is determined by a variational principle, see 0[7]. 



4.2 Renormalizability of the quark-gluon model in the presence of the Gribov- 
Zwanziger term 

In order to discuss the renormalizability of expression (|79p , one uses the Ward identities that have already 
been established in [31 1211 [5| [6] 17]. Moreover, it exhibits the following symmetry 

e^Etot = , (82) 
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where the operator is given by 

O\0,b f ( ,„ac 8 -be i , ,ac $ -be $ , jrac $ frac & 

+ U ^7f7be- ~ ^JffbA • (83) 

The Ward identity (|82j) expresses the invariance of the action (j79[) under a global U(A(N 2 — 1)) transfor- 
mation. This symmetry works exactly as the global ?7(4iV) symmetry associated with the spinor sector 
of the theory through the operator (|23|) ■ These global symmetries ensure in fact that no mixing terms be- 
tween the two set of BRST doublet fields, i.e. (<^f , y^f , wJJ 6 , co^ b ) and (C a , Oa,Vai Ki)i ar i se i* 1 tne allowed 
counterterm. All Ward identities of the Gribov-Zwanziger action remain valid in the present case. This 
is also the case of the identities ([28]) - ([38]) . Of course, the Slavnov- Taylor identity ([28]) is supplemented 
by suitable extra terms accounting for the BRST new doublets of the gluon sector, (|7ip and (|75p . 

5(E) - 5 (E„) + / A + + + . (84) 

Also, the ghost equation (|32p needs a little modification which, due to the presence of the Gribov- 
Zwanziger term, generalizes to 

g ^ r + ,/* / a + <^ + + e-^) , (85) 

while the classical breaking term (|34p remains unmodified. The previous algebraic analysis can be now 
repeated for the more general action (|79p . The final output is that the action (|79p remains renormalizable 
to all orders. 

5 Conclusion 

In this work we have considered a model that accounts for a modification of the infrared behavior of quark 
and gluon propagators in Yang-Mills theories. This is achieved through the introduction of suitable mass 
parameters which give rise to a soft breaking of the BRST symmetry, as outlined in |13| . 

Being soft, the breaking term can be neglected in the ultraviolet region, where the standard mass- 
less quark propagator is recovered as well as the notion of exact BRST invariance. Moreover, in the 
infrared region the quark propagator turns out to be deeply modified, as shown by expression ([2]). The 
physical reasoning behind the introduction of the soft BRST breaking and of the ensuing modification of 
the propagator relies on quark confinement and on the breaking of the chiral symmetry, both occurring 
in the non-perturbative infrared region. It is worth remarking that the quark propagator ([2]) is in fact in 
qualitative agreement with the fitting formulas employed in the numerical studies of the quark two-point 
function through lattice simulations in the Landau gauge [T71 [IB] . 
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The main result of the present article is the analysis of the renormalizability of the model, which we 
have shown to hold at any given finite orders of perturbation theory, by making use of the algebraic 
renormalization |19j . The inclusion of the Gribov-Zwanziger term which enables us to implement the 
restriction to the Gribov region Q has also been taken into account. Despite the presence of the soft 
BRST breaking term, the renormalizability of the model is guaranteed by the large set of Ward identities, 
eqs. ([28|) - ([38|) . which can be established. 

We expect that the mechanism of introducing non-perturbative infrared effects through the soft break- 
ing of the BRST symmetry [13] applies as well to other kinds of models, including supersymmetric and 
topological field theories. 
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